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Abstract. We present the type la rate measurement based on two EROS supernova search campaigns (in 1999 
and 2000). Sixteen supernovae identified as type la were discovered. The measurement of the detection efficiency, 
using a Monte Carlo simulation, provides the type la supernova explosion rate at a redshift ~ 0.13. The result is 
0.125±g;g|4± ; o l| /i 2 SNu where 1 SNu = 1 SN / 10 10 L| / century. This value is compatible with the previous 
EROS measurement iHardin et al.ll2000h . done with a much smaller sample, at a similar redshift. Comparison 
with other values at different redshifts suggests an evolution of the type la supernova rate. 
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1. Introduction 

Type la supernovae are believed to be thermonuclear ex- 
plosions of a white dwarf reaching the Chandrasekhar 
mass after accreting matter in a binary system (see e.g. 
iLiviol 1120011) ). They are the main mechanism to enrich the 
interstellar medium (ISM) with iron-peak elements (with 
~ 0.6 M© of nickel released per event on average). The 
knowledge of the number of such events per galaxy and 
per unit of time is crucial to understand the matter cycle 
and the ISM chemistry. From a cosmological point of view 
the measurement of the evolution of the explosion rate 
will put strong constraints on galaxy evolution. Measuring 
the evolution of the core-collapse supernovae (type II and 
Ib/c) explosion rate would give a new unbiased view on 
the star formation rate (SFR) history, independent of any 
peculiar tracer, because gravitational supernovae are di- 
rectly linked to massive star formation. Up to now, the 
type II super nova rate has been mea sured only in the lo- 
cal Universe l|Cannellaro et al.lll999h . On the other hand, 
the measurement of the evolution of the explosion rate of 
type la supernova would give strong constraints on the 
SFR history and on many other parameters such as the 
progenitor evolution (birth rate of binaries, time delay be- 
tween the white dwarf formation and the thermonuclear 
explosion...) or the nature of the progenitor itself. 

This paper is organized as follows: first the EROS 
search for supernovae is reviewed in Sect. Then the 
two supernova search campaigns at the origin of this work 
are presented in Sect. [31 Sections 0] to deal with the rate 
computation. Our results are presented in Sect. [Bland dis- 
cussed in Sect. ED 

2. The EROS search for supernovae 

The EROS experiment used a 1-meter telescope based 
at La Silla observatory, Chile, designed for a baryonic 
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scope beam light was split by a dichroic cube into two 
wide field cameras (1 deg 2 ). The cube sets the passband 
of both cameras, a "blue" one overlapping the V and R 
standard bands, another "red" one matching roughly the 
/ standard band 1 . Each camera was a mosaic of eight 2048 
x 2 048 pixels CCDs with a pixel size of 0.60" on the 
sky l|Bauer fc de Katl Il997t IPalancme-Delabrouille et, alJ 
1998). In the following, we will use only the blue cam- 
era. One "image" stands for one CCD, and one "field" is 
for the entire camera with its eight CCDs, i.e. ~ 1 deg 2 . 

A wide field camera is a good asset to perform super- 
nova searches as the number n of detected supernovae is 
proportional to fl ■ zf im where f2 is the solid angle surveyed 
and zi; m is the limiting redshift of the search - typically 

Send offprint requests to: G. Blanc, e-mail: 
blanc@pd . astro . it 

By "standard" we mean the Johnson-Cousins U BVRI pho- 
tometric syste m jjohnsonl Il965l : ICousinsI Il976l) as used by 
lLandold (Il992ft . 



EROS detected type la supernovae up to z = 0.3. EROS 
dedicated about 15 % of its time to supernova search be- 
tween 1997 and 2000. The main limitation has been to 
secure enough telescope time to perform the photometric 
and spectroscopic follow-up of discovered events. 

To search for supernovae EROS uses the CCD subtrac- 
tion technique which consists in observing the same fields 
three weeks to one month apart in order to detect tran- 
sient events and catch only supernovae near maximum 
light (a typical type la supernova reaches its maximum 
light in the optical wavelength roughly 15 to 20 days after 
the explosion). Both observations are done around new 
moon in order to limit the sky background. The two sets 
of images are then automatically subtracted after a spatial 
alignment, a flux alignment and a seeing convolution. 

Cuts are applied to eliminate known classes of variable 
objects such as asteroids, variable stars, and cosmic rays. 
The most important cut is the requirement that the tran- 
sient object be near an identified galaxy. More precisely, it 
is required to be within an ellipse centered on any galaxy 
with semi-major and minor axes of length 8 times the 
r.m.s. radius of the galaxy's flux distribution along each 
axis. This ellipse should contain essentially all supernovae 
associated with a galaxy since ~ 99 % of galactic light 
is contained within an ellipse of size 3 times the r.m.s. 
radius. To eliminate cosmic rays , each exposure consists 
in two subexposures: SExtractor ijBertin fc Arnout3l 996) 
software is run on both images to make a cosmic-ray cat- 
alog before adding the two images. 

In spite of these cuts, a human eye is still needed to 
check the remaining candidates and discard obvious arte- 
facts (residual cosmic rays, bad subtractions, asteroids, 
variable stars...). Ten to fifteen search fields are observed 
every night and the human scanner has to deal with about 
100 candidates the following day. The serious candidates 
are observed again on the following night for confirmation. 
Then a spectrum is taken to obtain the supernova type, 
phase and redshift. 

EROS supernova search fields lay near the celes- 
tial equator (declination between —4.5° and —12°) to 
be reached from both hemispheres and between 10 h 
and 14. 5 h in right ascension. They over lap some of the 
Las C ampanas Redshift Survey (LCRS) l|Shectman et alJ 
1996jj) fields which will facilitate the galaxy sample cali- 
bration. 

3. The 1999 and 2000 campaigns 

At the beginning of 1999, the Supernova Cosmology 
Project (SCP) undertook a large nearby supernova search 
involving eight supernova search groups, including EROS. 
The goal was to gather a large set of well sampled, CCD 
discovered, nearby [z < 0.15) type la supernovae in order 
to calibrate the distant events used to measure the cos- 
mological parameters. The whole collaboration obtained 
37 SNe out of which 19 were type la SNe near or before 
maximum light and were followed both spectroscopically 
and photometrically. EROS observed 428 square degrees of 
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SN 2000bw 
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20.5 
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SN 2000bx 
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Table 1. List of supernovae discovered by EROS in 1999 and 2000. From left to right: IAU name of each event, 
discovery date, discovery V magnitude, SNe celestial coordinates, redshift, type, phase (discovery time from maximum 
light) and IAU. 



sky in two steps and discovered 12 type la (cf. IAUC 7117 
and 7136) of which 7 have been followed by the collabo- 
ration. After obtaining reference images of search fields in 
January 1999 (01/12 to 01/29), the SN search was con- 
ducted between February 4 and February 27 (new moon 
was on February 16), and between March 09 and March 
27 (new moon was on March 17). 

One year later, EROS was involved with the European 
Supernova Cosmology Consortium (ESCC, involving 
French, British, Swedish and Spanish institutes) to search 
for intermediate redshift (z ~ 0.2 — 0.4) supernovae. Four 
type la supernovae were discovered (cf. IAUC 7406). One 
hundred and seventy square degrees were observed be- 
tween 2000, March 27 and April 9 (new moon on April 
4); reference images were taken between February 27 and 
March 15. The main difference between the two searches 
was the exposure time used: for the 1999 campaign it was 
set to 300 seconds, while for the 2000 search it was 600 
seconds allowing a deeper search by 19 % in redshift (as 
zi; m ex ToVp , where T cxp the exposure time). 

In both campaigns, spectra of each candidate were 
taken. The main characteristics of the discovered super- 
novae are summarized in Tabled 

4. Principle of the measurement 

Our supernovae search requires that a supernova be as- 
sociated with an identified host galaxy. We must there- 
fore make an assumption about how the supernova rate 



scales with galaxy luminosity, so as to correct the rate 
for supernovae in dim, unidentified galaxies. We choose to 
assume that the number of supernovae per galxay and 
per unit time is proportional to the red galaxy lumi- 
nosity, an assumption that receives so me empirical sup- 
port, though mostly in the blue band l)Tammannlll970i 
ICappellaro et allll993|) . Under this assumption, the ex- 
plosion rate TZsn is given by the ratio of the number of 
detected supernovae of type la A^sNe, to the number of 
galaxies to which the search is sensitive A/" ga i , weighted by 
their mean luminosity (£ ga i) and by the mean time in- 
terval (T) during which the supernovae can be detected: 

The time interval or control time T during which a given 
supernova is visible depends on the supernova detection 
efficiency e. Since the efficiency itself depends on many 
parameters, the observed control time T for a supernova 
of redshift z is given by an integral of e over these param- 
eters: 

/• r+oo 

T(z)= / / e(A; m(t, z; </>sNia)) dt d^sNia, (2) 

J 0SNIa J — OQ 

where A is the time interval between the search images 
and the reference images; m(t, z; </>sNia) is the supernova 
apparent magnitude (t is the time from maximum light 
or phase; z is the redshift). The magnitude m depends 
also indirectly on the light curve distribution of type la 
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supernova (see S I7.3[) . that we note 0sNia- As we know the 
galaxy sample in the search images (see §[EJ, but not their 
redshifts, the restframe explosion rate can be rewritten as: 



A/sNe 



E ff Jo + °° P( z \ m g) ■ c g( m 9i z ) ■ r ( z )l (l + z)-dz' 



(3) 

where g is a subscript running on all galaxies of all the 
search images; p(z\m g ) is the probability density of the 
redshift z of a galaxy knowing its apparent magnitude 
m g (see § I7.2fl : C g is the galaxy luminosity (which de- 
pends on its redshift and apparent magnitude, as well 
as on a cosmological model: we use a standard model, 
(flfcfoj ^A ) = (0.3, 0.7); for the Hubble constant, we use 
the notation H a — 70 ft.70 km s _1 Mpc -1 wherever 
needed). Different steps are needed to calculate the de- 
nominator of eq. © . The first one is to establish a well- 
defined sample of galaxies in the search images. This is 
then used to measure the detection efficiency. Each of 
these steps is detailed in the following sections. 

In the following the EROS instrumental magnitude 
will be denoted: 



#cros = -2.5 • log 



Fadu 



(4) 



cxp 



where Tadu is the flux as measured in CCD images 
in digital units; T C xp is the exposure time in seconds. 
Throughout this paper the magnitudes in EROS images 
are computed using an adaptive apert ure photometry as 
provi ded by the SExtractor software i Bertin fc Arnoutl 
which is accurate enough for our purpose. 
Photometry for stellar objects was done using observa- 
tions of Landolt stars, as detailed in Sect. 17.41 Galactic 
magnitudes needed to derive the redshift distribution and 
to normalize the SN rate was calibrated using LCRS R 
magnitudes. 



5. The galaxy sample 

We use the SExtractor software l|Bertin fc Arnoutsl ll996) 
to detect objects on the images. Objects are classified as 
stars and galaxies using the SExtractor CSTAR parameter 
which quantify the "stellarness" of an object. For stars 
CSTAR ->• 1, and for galaxies CSTAR -> 0. The critical 
input in order to obtain a reliable value of this parameter 
is the image seeing. SExtractor is then run twice on the 
images, first to obtain the seeing, and a second time to 
compute CSTAR. Figure [T] shows a scatter plot of CSTAR 
versus the galaxy magnitude. We see that the galaxy star 
separati on is clear for R < 19 . The galaxies of the LCRS 
catalog l| Sheet man et al . 1996) shown in the plot allow one 
to verify the efficiency of the classification for R < 18. 

To further check the reliability of the galaxy sample 
we make a "galaxy count" that we compare to published 
ones. As illustrated in Fig. [5] we are able to define ac- 
curately the completeness limit of our survey as the point 
where our data turn down relative to the published counts. 
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Fig. 1. EROS magnitude distribution vs. value of CSTAR. 
Objects with CSTAR around 0.5 are faint and thus dif- 
ficult to classify. Squares stand for LCRS galaxies. The 
dark line represents the cut © for the 1999 search. 



At this point, stars and galaxies are quite hard to distin- 
guish: most of the objects still have a CSTAR ~ 0.5 (see 
Fig. n]for B croa ~ —2.6). To take this into account, we set 
the galaxy magnitude cut one magnitude below the com- 
pleteness limit in order to prevent contamination of our 
galaxy sample by stars. Then: 



for 1999 



CSTAR < 0.55, 

B cms < -3.5; {i.e. R<19), 



(5) 



for 2000 



(6) 



CSTAR < 0.55, 
Scros < -3-1; (ie. R < 19.3). 

The effect of this cut is shown in Fig. [T] We used the 
LCRS galaxy catalog ijShectman et al.lll996h to check it, 
as illustrated in Fig. ^ This cut removes events in faint 
galaxies which are visually indistinguishable from variable 

stars. 

According to the LCRS luminosity function ijLin et alJ 

Il996|) . our cut requiring R < 19 means that at our mean 
redshift of 0.13 (see Fig. fTTft . our supernova search is sen- 
sitive to supernovae in galaxies that generate about 85 % 
of the total stellar luminosity. We are not sensitive to su- 
pernovae in the multitude of low-luminosity galaxies but 
these galaxies generate little total light. If the supernova 
rate is proportional to the luminosity, they also generate 
few supernovae. 

6. The detection efficiency 

In order to measure the detection efficiency, we super- 
impose simulated supernovae on the galaxies. On each 
search image we choose a bright - but not saturated - 
star (with a signal to noise ratio typically above 50) as a 
template to be "copied" onto the potential host galaxy. 
Among all such stars found in the image, we take the 
one nearest to the galaxy on which we wish to add a 
supernova. We move only the pixels included in a circle 
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(a) Galaxy count on EROS fields done on 180 images 
from the 1999 search. The completeness limit is for 
CSTAR < 0.55 and S OI . OS < -2.6. 



(b) Galaxy count on EROS fields done on 29 images 
from the 2000 search. The completeness limit is for 
CSTAR < 0.55 and Seros < -2.2. 



Fig. 2. Galaxy counts in EROS supernova fields for various cuts on the CSTAR parameter. The nominal value for 
this cut is CSTAR < 0.55. Beyond, faint stars pollute th e galaxy sample. Below, fai nt galaxies are missed. The curve 
stands for published data such as the count in R band bv lBertin &: Dennefeldl <|1997l) . The matching of the curve with 
the histogram confirms our galaxy magnitude calibration, eq. H13|) . The main difference between the 1999 and the 2000 
searches was the exposure time, set to 300 seconds in 1999 and 600 seconds in 2000. 



with a radius given by the point where the star flux cal- 
culated by using a gaussian is equal to half a standard 
deviation of the sky background. Each pixel is rescaled to 
the desired flux (simulated supernova flux is distributed 
uniformly between and 15 000 ADU - the latter corre- 
sponding roughly to V ~ 18.5 for 1999 and to V ~ 19.2 
for 2000). Simulated supernovae are placed at random po- 
sitions within the isophotal limits of the galaxy, i. e. inside 
an ellipse centered on the galaxy and that has semi-major 
and minor axes of length 3 times the r.m.s. of the galaxy's 
flux along the axis distribution; this ellipse contains ~99 % 
of the galaxy flux. The positions are chosen following a 
probability proportional to the local surface brightness. 
They are placed on the galaxies of the images before the 
search software is run. A simulated supernova is detected 
by the search software if its signal-to-noise ratio is greater 
than 5 and if it is situated within a radius of 1.5 pixels 
of the simulated one, according to Fig. |3] The resulting 
efficiency as a function of the simulated SN flux is the ra- 
tio between detected supernovae and simulated ones. The 
result is fitted by an analytic function in order to smooth 
the histogram. The efficiency used in eq. © is computed 
for each search image. A "global" efficiency computed for 
each of the two campaigns as a function of simulated flux 
is shown in Fig. 0] with the resulting fits. Figure shows 
this "global" efficiency as a function the simulated SN 
distances to their hosts. We remark that we observe no 
significant efficiency loss in the center of galaxies. This 
improvement comes from using CCD detectors instead of 
photograp hic plates (where th e core of galaxies saturates 
- see e.g. iHowell et al.l l|200fl(l ). Moreover the efficiency 
is roughly independent of the position within the central 



~ 8" of the galaxies. Beyond this radius, corresponding 
to the isophotal limit of most galaxies in our sample, few 
supernovae are simulated leading to the large statistical 
fluctuations in Fig. 03 Near the center of galaxies, the effi- 
ciency is about 0.9, limited mostly by "geometrical" losses 
near CCD edges and around dead pixels. 




Entries 


118316 


Mean 


0.5808 


RMS 


0.33 



3 4 5 

Distance simulated-detected (pixels) 

Fig. 3. Histogram of the distance between the simulated 
supernovae and the detected ones. The cut we used to 
identify the detected population is set at a radius of 1.5 
pixels from the simulated supernova as shown by the dark 
line. 
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Fig. 4. Mean efficiency computed for each of the two 
searches as a function of the simulated flux, in ADU s _1 
or V magnitude (according to eq. I|12jl ). for both the 1999 
search (continuous line) and the one of 2000 (dashed line) . 
The curves are the fits. With its exposure time twice the 
one of 1999, the 2000 search went about 0.2 mag deeper. 
The two vertical lines represent the limiting magnitudes 
defined as the magnitude at half the maximum efficiency: 
5.59 ADU s" 1 for 1999 and 5.09 ADU s" 1 for 2000. The 
ratio between the two limits is 1.10, to be compared with 
1.41 expected for an exposure time ratio of 2. The plateaux 
reached at high flux are respectively 0.94 and 0.91 for 1999 
and 2000. 



ffi 0.8 



0.6 



0.4 



0.2 



— 1999 
■ ■■ 2000 



domly drawn among a selected set of SNe (see Sect. 17. 
a redshift is drawn according to the probability density 
p(z\m g ) as defined in Sect. 17.21 the selected light curve is 
adjusted according to this redshift. A phase is drawn uni- 
formly, which provides a standard magnitude. After trans- 
formation of standard V and R magnitudes into ADU (see 
Sect. I7.4|l . the corresponding efficiency is read from Fig.QJ 
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Fig. 6. Examples of redshift distributions of galaxies 
knowing their apparent magnitude, as given by (J7J- The 
solid lines show the predicted distribution for three ap- 
parent magnitudes m g and two sets of cosmological pa- 
rameters. The histograms show the distribution of LCRS 
redshifts for galaxies within 0.1 mag of the given magni- 
tude. 



SN - Host galaxy distance (arcsec) 

Fig. 5. Mean efficiency as a function of the distance of sim- 
ulated supernovae to their host galaxy peaks, computed 
for each of the two searches. Since there is a loss of effi- 
ciency at low fluxes, this plot is done for fluxes above 20 
ADU s" 1 for 1999 and above 13 ADU s" 1 for 2000. 

7. The integral computation 

7.1. The algorithm 

To compute the integral in the denominator of eq. 
we have to integrate the detection efficiency over various 
parameters, such as the type la supernova light curve dis- 
tribution, the supernova redshift and its phase at the time 
of discovery (time from maximum light). In order to do 
that, we use a Monte-Carlo method: a light curve is ran- 



7.2. Redshift distribution of galaxies 

We do not know the redshifts of all the galaxies in the 
images. But all we need is the probability distribution p 
for the redshift z given the galaxy's apparent magnitude 
m g : 

, . . dN dV c . . 

where N is the number of galaxies having this redshift 
and apparent magnitude, Q is the solid angle surveyed, 
V c is the comoving volume and <p g the galaxy luminosity 
function {M g being their absolute magnitude). To com- 
pute this, we u se the luminosity function of the LCRS 
l(Lin et a"flll99f)l) . Figure El gives examples of such a distri- 
bution. 
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Fig. 7. Absolute V magnitude at maximum light distribu- 
tion for the whole published set (51 SNe, white histogram, 
upper box) and the reference sample of 17 supernovae we 
selected (blue histogram, lower box). If we remove the 
lowest luminosity event (SN 1996ai) from this reference 
sample, its mean is shifted from —18.55 to —18.75. 




Fig. 8. Distribution of the light curve shape parameter 
Amf 5 for the whole published set (51 SNe, white his- 
togram, upper box) and the reference sample of 17 su- 
pernovae we selected (blue histogram, lower box). 



7.3. Type la supernova light curve distribution 

We need to integrate the efficiency on a large and repre- 
sentative sample of supernova light curves, so as to repro- 
duce the variety of such objects. The light curve distribu- 
tion of type la supernova can be parametrized, at least at 
first order, using two parameters: the luminosity at max- 
imum light, and a light curve shape parameter. The light 
curve shape ca n be quantified using the Amf 5 parameter 
(Phillrpl llQQl . that is the magnitude difference between 
the light at maximum and 15 days later on the B light 
curve. 

Out of 51 published type la supernovae o bserved light 
curves ijHamuv et al.lll996tlRiess et al.lll999|h we selected 
a reference sample of seventeen objects that have a good 
sampling in both V and R bands, the standard bands used 
to calibrate the EROS band. Each light curve is then fitted 



with an analytic template l|Contardo et all ioOO) , which 
allows us to reach any point of the light curve without 
interpolation; the selected objects of the reference sample 
have a very good fit. The resulting distribution of absolute 
magnitude at maximum light (the luminosity function) is 
shown in Fig.0 The reference sample we have used has a 
lower mean luminosity than the full sample. This is mainly 
due to one very subluminous event, SN 1996ai (which may 
be ver y extinguished by i ts host galaxy - Ay ~ 4 accord- 
ing to lTonrv et ail HqQ^)). The Amf 5 distribution of the 
reference sample is shown in Fig. El it is similar to that of 
the full sample. Our supernova rate will be presented in 
a form such that it can be revised if future measurements 
give different values for the mean luminosities and Amf 5 . 

7.4. From standard magnitudes to ADUs 

As the detection efficiency is measured as a function of su- 
pernova flux in the EROS band, and as the published light 
curve sample we will use to reproduce the supernova light 
curve distribution is expressed in standard magnitudes, 
we need to translate one system into the other using a 
calibration relation such as: 



V + a-{V-R)+P + 8-(£-l), 



(8) 



between the measured flux on the CCD images (B CIOS is de- 
fined by eq. 10} ) and the observed magnitudes through the 
standard bands: we choose V and R which are the closest 
to the EROS band; a is the color term; f3 is the so-called 
zero point; 6 • (£ — 1) (£ is the airmass) is the atmospheric 
absorption: the atmospheric absorption at meridian (i.e. 
at ~ 20° from the zenith) is included in the zero point. 
V and R are related to the intrinsic magnitudes of the 
supernova V sn and R sn after correcting for the various 
extinctions: 



V = V sn + A, 



a extra 



R = R sn + A R + A°j 



(9) 



(10) 



where Ay and An are the Galactic absorption; Ay ctra and 
^cxtra are extragalactic absorption which could be ei- 
ther an intergalactic extinction or a host galaxy extinction 
or both. 

Moreover, we use light curves of observed nearby su- 
pernovae. These light curves are transformed to any red- 
shift; the relationship between the published light curve 
m F (t\,z{) in a filter F and the "new" redshifted light 
curve m F (t 2 , z 2 ) 2 is given by 



m F (t 2 , z 2 )-m F (t 1 , zi) = 5 log 



V{zi) 



K F (t 2 ,z 2 ) - K F {t\,Z\), 



(11) 



2 The supernova phase t\ at z\ is related to t-x as t\ / (1+Zi ) = 
t 2 /(l + z 2 ) 
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where T> is the luminosity distance and if is a term which 
takes into account the fact that we observe objects whose 
flux is sliding according to their redshift through a fixed 
passband in wavelength space. This term is known as the 
"if -correction" (see Sect. 17.4. 5f) . Note that this relation 
does not depend on the Hubble constant. The following 
sub-sections deal with the different terms described above. 



7.4.1. Can we calibrate supernovae using stars? 
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Fig. 9. Evolution with the supernova phase of the error 
we made when using stars to calibrate supernovae. What 
is actually plotted is the difference between zerojpoints 
computed for SNe la (B C tos — V) S ne and for stars (B CTOS — 
y) s tars- Each dot stands for a star whose color (V — R) is 
equal to the (V — R) color of a SNIa at that phase. The 
signal difference is never larger than 8 %. The colors are 
computed from template spectra. 

As SN spectra are very different from star spectra, 
showing broad blended features evolving along the SN 
phase, we may wonder to what extent we can use stars 
to calibrate SN fluxes. In order to quantify this effect, we 
use synthetic colors computed with template spectra for 
each phase of the SN. Figure shows the error we make 
when calibrating type la supernova fluxes using stars, as 
a function of the supernova phase. Assuming the color is 
the same for SNe and stars (asNc ~ a sta rs), the plotted 
quantity (# e ros - V) snQ - {B CIOS - V) s tais reflects directly 
the zero point differences, PsNe — /3stars- The error is at 
most 8 % around —10 days and between +10 and +20 
days. We note that for standard filters, there would be no 
such error. 



7.4.2. Zero point and color term for the stars 

Zero points for the EROS obs ervations are computed from 
standard stars l|Landoltlll992|) or secondary standard stars 
around some EROS super novae which were fol lowed-up 
and analyzed by the SCP l|Regnault et al.ll200lh . As the 
observations were always done close to the meridian, the 



atmospheric absorption correction is included in the zero 
point. The evolution of the zero point as a function of time 
is plotted in Fig. ED Within a 5 % accuracy, it is the same 
for the 1999 and the 2000 searches. The color term is com- 
puted using the SMC and LMC OGLE calibrated catalog 
ilUdalski et afll2000l) in order to have a larger lever arm: 
the Landolt catalog has the drawback to have mainly red 
(87 % have B — V > 0.3) stars. This can be overcome with 
the OGLE catalog of SMC and LMC stars, that EROS is 
surveying too in its microlensing search program. The re- 
sult is 

Soros = V - 0.53 ■ (V - R) - 22.7. (12) 

The zero point above is obtained for one CCD of the 
EROS mosaic, the seven other CCDs (that have different 
gains) are calibrated relative to this one, using the sky 
background and various photometric catalogs. The rela- 
tive variation from CCD to CCD is of order of 10 %. 

7.4.3. Zero point for the galaxies 

To calibrate the galaxy magnitudes we use the LCRS cat- 
alog which prov ides a single isophotal i?i cr s magnitude 
(|Lin et alJll996^ . We find: 

Soros = idlers _ 22.7. (13) 

iLin et alJ l|l99fj) estimate that R = Ri CIS - 0.25. The dis- 
tribution of Soros — -Ricrs for the LCRS galaxies has a r.m.s. 
width less than 0.2 mag. This width is sufficiently narrow 
that it was not necessary to refine eq. (|13fl with a color 
term. The reason that little precision is necessary is that 
the galaxy magnitudes are only used, first, to derive the 
redshift probability distribution for each galaxy and, sec- 
ond, to calculate the total galactic luminosity of the sam- 
ple to normalize the supernova rate. In both cases, one 
effectively averages over all galaxies in the sample so the 
relation lfLl|) is sufficiently accurate. 

7.4.4. Extinction 

Observations are always performed close to the meridian, 
so the resulting atmospheric absorption can be considered 
constant and included in the zero point. 

The Galactic extinction is corrected field by field us- 
ing the ISchlegel et ail l)l998l) reddening map, assuming a 
standard Galactic extinction law. The mean absorption is 
A v ~ 0.12 for the EROS fields. 

We do not correct for the supernova host galaxy extinc- 
tion (or intergala ctic extinction) since the templat e light 
curves from lHamuv et alJ (^996) and iRiess et all l|l999t) 
do not correct for it either. We thus implicitly assume 
that the absorption in these two samples is representative 
of supernovae in general. 

7.4.5. if -corrections 

Since the supernovae are observed at different redshifts, 
their spectra will move across the fixed passband used. 
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(a) Zoom on the 1999 supernova search: the mean 
zero point is about 22.7 during the supernova search 
(see first box for the February campaign, and second 
box for the March campaign). 



(b) Zoom on the 2000 supernova search: the mean 
zero point is 5 % lower than the one from 1999 
searches (see box). 



Fig. 10. Evolution of the "blue" zero point as a function of time. Circles stand for the calibration using the Landolt 
catalog; black filled circles are for secondary calibrating stars from the supernova calibration of SCP 1999 set; crosses 
are absorption coefficients from EROS analysis of LMC and SMC stars, shifted to match the Landolt points. The 
February 1999 search has been done between 02/04 and 02/27; the March 1999 search was performed between 03/09 
and 03/27; The April 2000 search took place between 03/27 and 04/09. For these three campaigns the corresponding 
zero point is 22.7 within a 5 % range. The total uncertainty on this measurement is, from the scatter of the points, 
about ±7 %. 



In a given passband, the sampled part of the spectrum 
thus depends on redshift and this influences the measured 
flux. This is taken into account by the "K -correction" . 
The correction in the band F is given by 



K F (t,z) = -2.5 -log 



J+°°d\ F{A )I(t,\ ) 



(14) 

where A Q is the observed wavelength, F is the filter trans- 
mission function and / is the normalized time-dependent 
spectrum, which may depend on time t - like for super- 
novae. Note that Kp = when F — 1, i.e. when the 
passband covers the whole spectrum. 

For type la su pernovae they are c alculated using a 
spectral template 3 I Nugent et alJ l2002) and the shape of 
the standard V and R filters (which are the bands used 
to calibrate the EROS band, according to eq. (JUJ). For 
galaxies we use the approximation K r(z) = 2.5 -log(l + z) 
which is good enough up to z ~ 0. 2 |Poggiantilll997|) . 



3 If we use the improved spectral template bv lNobili et alJ 
(2003) the resulting value of the rate is increased by ~ 1 % 
which is negligible. 



8. Results 

To check the validity of the Monte Carlo simulation we 
compare some simulated distributions to the observed 
ones. Quantitative comparisons are performed using a 
Kolmogorov-Smirnov (KS) test. Figures ITT1 to IT31 give ob- 
served and simulated distributions for various parameters, 
namely the redshift, the magnitude at detection and the 
host magnitude. The similarity between the observed and 
the simulated distributions gives confidence in the relia- 
bility of the integral computation. Considering the limited 
statistics and the KS probabilities (85 % for the redshift 
- Fig. ^] 62 % for the discovery magnitude - Fig. 
23 % for the host galaxy magnitude - Fig. fH3|) give no ev- 
idence for systematic differences between theoretical and 
observed distributions. 

The distribution of the distance between the super- 
nova and the host center is shown in Fig. The Monte 
Carlo distribution was generated assuming that the super- 
nova rate is proportional to the local surface brightness. 
The observed distribution is more concentrated toward 
the center than the surface brightness with the KS test 
giving a 3 % probability for compatibility. This suggests 
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Fig. 11. Observed (red, filled - upper box) and simulated 
(blue, hatched - lower box) redshift distribution. The two 
empty bins are SNe for which the host galaxy magnitude is 
fainter than cuts JHJ or JJjJl. A Kolmogorov-Smirnov test 
gives a 85 % probability of adequacy between the two 
distributions. 
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Fig. 13. Observed (red, filled - upper box) and simu- 
lated (blue, hatched - lower box) host galaxy R magni- 
tude distribution. The two empty bins are SNe for which 
the host galaxy magnitude is fainter than cuts © or ©. 
A Kolmogorov-Smirnov test gives a 23 % probability of 
adequacy between the two distributions. 
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Fig. 12. Observed (red, filled - upper box) and simu- 
lated (blue, hatched - lower box) discovery V magnitude 
distribution. The two empty bins are SNe for which the 
host galaxy magnitude is fainter than cuts or |JBJ. A 
Kolmogorov-Smirnov test gives a 62 % probability of ad- 
equacy between the two distributions. 



supernovae do not have the same distribution as the light 
though the limited statistics prevent any firm conclusion. 
At any rate, except far from galactic centers, our detection 
efficiency is nearly independent of the distance to the host 
center so the precise distribution of the distance has little 
effect on the global efficiency. All type of host galaxies are 
included in the simulated distribution. The observed dis- 
tribution confirms that we are able to detect supernovae 
in the center of their host galaxies. 

These checks being done, we are able to compute the 
rate, as the ratio of the number of detected supernovae 
to the integral value calculated as above. From the 1999 
search we have 12 type la supernovae. But cut (JSJ) on the 
galaxy host magnitude discards two events (SN 1999bl 
and SN 1999bn). In 2000, four type la supernovae were 
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Fig. 14. Observed (red, filled - upper box) and simulated 
(blue, hatched - lower box) distance to host distribution. 
The two empty bins are SNe for which the host galaxy 
magnitude is fainter than cuts JSJ or A Kolmogorov- 
Smirnov test gives a 3 % probability of adequacy between 
the two distributions. The simulated distribution assumes 
that the type la supernova rate is proportional to the 
galaxy surface brightness. 



detected. Only 14 events remain after the cut. We can 
rewrite eq. (j2J as 



— 



(15) 



where A/tot is the total number of events surviving the 
cuts, and Itot is the corresponding value of the inte- 
gral. For the 1999 search, the integral is I" = 0.89 x 



10 14 hZj Lg yr, while it is X 00 = 0.50 x 10 1 



h 7 Q Lg yr 

for the 2000 search. The combined value of the integral 
for the two searches is the sum X tot = X 99 + T m = 
1.39 x 10 14 /iy 2 Lq yr. Assuming the number of ob- 
served supernovae follows a Poisson law, with a confi- 
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dence level of 68.3 %, one has Mot = 14±l;™. This 
gives Kf a = O.lOOtg-.ga? ^70 SNu fl , where 1 SNu fi is 
1 supernova/ 10 10 Lq /century. The errors are only sta- 
tistical at this point. The mean redshift for this value is 
given by the mean of the simulated redshift distribution 
(Fig. El, that is (z) = 0.128. The calculated r.m.s. width 
of this redshift distribution is 0.056. 

8.1. Systematic uncertainties 

We have identified four possible sources of systematics: 
the calibration for both supernovae and galaxies, the cos- 
mological model and the assumed distributions of super- 
nova luminosities and light-curve shapes. For each of these 
parameters, various simulations have been performed to 
quantify the corresponding impact on the rate. 

8.1.1. Calibration 

We estimate the uncertainty on the zero point of eq. (|12|1 
to be of 0.07 mag (see Fig. An increase (resp. de- 
crease) in the zero point of 7 % decreases (resp. increases) 
the rate by 6 %. 

Considering the galaxy calibration, eq. (| 1 31) . we esti- 
mate an uncertainty on the zero point of about 0.1 mag 
(as combination of statistical error in the ensemble zero- 
point error along with a possible systematic calibration 
error). An increase (resp. decrease) of the zero point by 
10 % increases the rate by 7 % (resp. decreases the rate 
by 6 %). The net effect is lower than the 10 % due to 
luminosity alone, because the calibration for the redshift 
distribution acts in the opposite way 4 . 

8.1.2. Cosmology 

The calculations and simulations have been done with the 
(0^f o , fl\ ) = (0.3, 0.7) cosmological model. We checked 
that using the (^a/ , ^a ) = (1, 0) model lowers the rate 
by about 1 %. Only the probability density p(z\m g ), the 
galaxy luminosities and the supernovae magnitude (see 
eq. fTT|0 . through the luminosity distance, depend on the 
cosmological model. But while the luminosity distance in- 
creases from (1, 0) model to (0.3, 0.7) model, on the con- 
trary, the p distribution is shifted toward the lower red- 
shifts for a given apparent magnitude. The two effects can- 
cel in this redshift range, hence the very small dependence 
of the cosmological parameters on the rate. 

8.1.3. Supernova diversity 

To take into account the fact that the distribution of su- 
pernovae luminosities and light-curve shapes is currently 

4 This can be understood in the following way: if the zero 
point is increased, the corresponding R magnitude increases 
for a given EROS flux. Then the redshift distribution is shifted 
toward higher redshifts. This change increases the value of the 
integral and thus decreases the rate. 



uncertain due to a lack of statistics, we parametrized the 
rate with the two main variables describing the type la va- 
riety distribution, the mean absolute magnitude at maxi- 
mum ( A/y +5 log hm) , and the light curve shape param- 
eter (Amf 5 ) (|Hamuv et al.lll99rj() for the 17 light curves 
we have used to derive the rate (see Fig. [7| and [SJ) . Then 
we compute the following empirical relation: 

( Kia \ = / (Amf 5 ) \° 65 
V 0.100 hj Q SNufJ V L138 / ( 16 ) 
x [1 + 0.78 • ((M Vmax + 51og/i 70 ) + 18.55)] , 

for deviations up to ~ 20 % for (My +51og/i7o) and up 
to 10 % for (Amf s ). If the whole type la variety distribu- 
tion is shifted, either in luminosity or in light curve shape 
parameter, as statistics increase with future experiments, 
our rate value will be able to be adjusted accordingly to 
eq. (Q3J). 

The uncertainty on the mean of the peak luminosity 
distribution (Fig. [TJ) of ~ 0.96/V17 = 0.23 mag translates 
into an uncertainty on the rate of 0.78 x 0.23 = 0.18. 
The uncertainty on the mean value of Amf 5 (Fig. |HJ|, of 
0.257/VT7 = 0.062, translates into an uncertainty on the 
rate of 4 %. Then, the main source of uncertainty on the 
rate comes from the mean value of the peak luminosity 
distribution of supernovae. 

If we remove the lowest luminosity supernova 
(SN 1996ai - see Fig. 0) from our reference sample, the 
rate is decreased by 16 %. 

8.1.4. Supernovae in dim galaxies 

Our search strategy requires an association to a host 
galaxy. This prevents detecting supernovae appearing in 
low brightness galaxies whose magnitude is beyond the 
detection threshold. For example, with such a cut, the 
type la SN 1999bw supernov a, the host of which i s a very 
low brightness dwarf galaxy ijStrolger et alJl2002|) . would 
not have been detected by EROS (this supernova had a 
B magnitude of 16.9 at maximum light, while its host has 
a B magnitude of 24.2). iGal- Yam et all l|2003h also found 
two hostless type la supernovae in galaxy clusters. They 
argue that these events can be due to a putative inter- 
galactic star population. 

Since we assume that the number of supernovae is pro- 
portional to the host galaxy luminosity, the derived rate is 
independent of the galaxy luminosity cut (eq. (JSJ) and ©) 
(apart from statistical fluctuation). Then such an eventu- 
ality is included in our rate derivation as long as our work- 
ing hypothesis is valid. Moreover we remark that these 
events are very rare, most probably representing a small 
fraction of events in the field. 

8.2. The rate 

Table El summarizes the error budget as discussed above. 
Taking it into account gives for the rate, the value: 

= 0.100t°;°r 7 ±° ;^ h% SNu«, (17) 
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Calibration: 



SNe 
galaxies 



SNe light curve distribution: (Amf 5 ) 

(M w ) 



±6 % 

+7 % 
-6 % 

±4 % 

±18 % 



Total (quadratic sum) ±21 % 

Table 2. Summary of identified systematic errors. The 
resulting total systematic uncertainty can be compared to 
the statistical uncertainty: _ 23 y 



where the first quoted error is statistical and the second 
systematic. This number comes naturally in the R band. 
However, the SN rate is traditionally expressed per lumi- 
nosity unit in the B band. We can convert this measure- 
ment into the R band as: 



tR 



= lo -Q.4.[(B-fl) -(B-Ji)g 8 i] i 



(18) 



■'al 



Using (B — R)q = 1.06 and ( B — R) va ] = 1.3 a s the m ean 
of the whole galaxy sample of iGrogin fc Gellerl (|1999f> : we 
take ±0.1 mag as a possible systematic between our galaxy 
sample and theirs. Then, we have 



tR 



T B 



1.25 ±0.11, 



which gives: 



V B - i25+ 044 + 028 h 2 SNu 

'Ma — u - izo -0. 034-0. 028 "70 oi, "U, 



(19) 



(20) 



where the uncertainty which appears in the right hand side 
of eq. (|19|l has been added quadratically to the systematic 
error. 

To obtain the rate per comoving volume unit, we mul- 
tiply the value in eq. (|20|l by the mean gala xy luminos- 
ity de nsity. The 2dF Galaxy Redshift Survey l|Cross et alJ 
200l| provides a value at a mean redshift of 0.1, in 
the B band and for a (1, 0) cosmological model, in 
agreeme nt with the one obtained by the ESO Slice 
Project l|Zucca et al.lll997(l . We translate the 2dF value 
to a (0.3, 0.7) model 5 ( j (0 ' 3 ' °- 7) /.7 (1 ' 0) | 2=0 : = 0.91): 
j B = (1-59 ± 0.13) x 10 8 h 70 L Mpc -3 . The super- 
nova rate Il2()[l can then be expressed as: 



(1.99l^r 4 4 7) x 10- 5 h 70 yr^Mpc- 3 . 



(21) 



Here, the uncertainty on the luminosity density function 
has been added quadratically to the systematic error bar. 

9. Discussion 

9.1. Comparison with other measurements 

Table |3| shows other published results for the type la 
supernova rate at various redshifts. Our measurement 

5 Note that j oc D 2 /V where T> is the luminosity distance 
and V = dV I ( dzdQ) is the comoving volume element (see e.g. 
ICarroll et all (Il992h . eq. (26)). 



compares well with other measur e ments in the nearby 
universe llCappellaro et alJ Il999l lHardin et all l200ol 
iMadewick et al.ll2003|) . It is slightly lower but compatible 
within error bars. We note that our SN search strategy 
and our supernova samp le a re very different from th ose of 
ICappellaro et alJ (|l999ji and lMadewick et all l|2003|) . The 
former used a sample of nearby supernovae discovered ei- 
ther by eye or on photographic plates which introduced 
large systematic errors in the rate. The latter used a new 
search technique based on galaxy spectra: supernovae are 
discovered in subtracting 116 000 galaxy spectra from the 
Sloan Digital Sky Survey by a eigenbasis of 20 "unpol- 
luted" galaxy spectra. The method is very promising, es- 
pecially in deriving the rate as a function of host galaxy 
properties. 

On the other hand our rate is lower than the distant 
IPam et all 0002) value, by over one standard deviation, 
either for the rate as a function of galaxy luminosity or as 
a function of comoving volume. 



9.2. Evolution? 

Our resu lt combined with the SCP measurement at higher 
redshift <|Pain et alJ l2002) suggests models with an evolu- 
tion of the rate, either in SNu or per unit of comoving vol- 
ume, between z — 0.13 and z — 0.55. We can parametrize 
such evolution as a power-law: 7?.i a (z) = 7?-i a (0) • (1 + z) a 
where a is a r estframe evolution index. Using the present 
result and the lPain et al.l l)2002h value per unit of comov- 
ing volume, we obtain a v — 3.1 



+ 1.6 



,1 4 . The same parameter 
computed for the rate expressed in SNu is ai — 2.6l" 1 ; 5 i . 
The difference, a v — ai = 0.5 ± 0.1, simply refle cts the dif- 
ferent luminosity densities adopted by us and IPain et alJ 
ll2002l) . equivalent to an evolving luminosity density with 
a fs ° = 0.5±0.1. In fact, the luminosity density most likely 
evolve s a bit faster. Combining the d ata of iLillv et alJ 
l|l996|) with that of ICross et all l|200lj) and imposing a 
(0.3, 0.7) cosmological model, the evolution of the lu- 
minosity density corresponds to a^ lly = 0.73 ± 0.51. If 
IPain et alJ §002) had adopted a luminosity density more 
in line with this value of aj B then we would have de- 
duced ai = 2.4+ } £ by combining our data with that of 
IPain et al.ll|2002l) . This number is still significantly greater 
than zero. 

The value of a„ = 3.1 + -. ' ^ deri ved by combining our 
data with that of lPain et al. 1 2002) is consiste nt with the 
value a„ = 0.8 ± 1.6 derived using only the IPain et alJ 
(2002) data. Nevertheless, fitting all available rate mea- 
surements either in SNu or per comoving volume unit, pro- 
vides a lower evolution index (see Fig. ll5f) : a v = 1.75±0.84 
and ai = 1.21 ±0.80. 

Models s uch a s those by iMadau et alJ l|l998l) or by 
ISadat et aD lll998l) can fit the observations (Fig. ^J. In 
the lMadau et alJ l)l998l) models, the evolution depends not 
only on the time delay between the white dwarf formation 
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Table 3. Comparison with other published restframe type la supernova explosion rate measurements. These rates are 
given together with the mean redshift of the observed SNe and the cosmological model assumed (especially for distant 
values). The fifth column gives the number of supernovae from which the rate is computed. Notes: a) the value per 
comoving volume unit is derived from the SNu value using the 2dF lumino sity density (the erro r on the luminosity 
density is added quadratically to the systematic uncertainty if distinct); b) lHardin et all l)2000|) computed the rate 
using a (Qu > ^A Q ) = (0.3, .0) cosmologica l mode l. Nevertheless no important difference from the (0.3, 0.7) model is 
expected (less than 10%); c) ICal-Yam et all 1 2002ft have measured the rate in galaxy clusters, with could be different 
from the rate in the field; they do not give a value per unit of comoving volume. 
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mation scenario, especially at high redshift. In these mod- 
els the supernova rate follows roughly a (l + z)" 11 evolution 
with a v > 1.8 fr = 0.3), q „ ~ 2.4 (r = 3) up to a v ~ 2.8 
(t = 1). The ISadat et all l|l998|) models depend on the 
choice of star formation rate evolution model and on dust 
extinction. They predict a v ~ 1.4 for "Ml" model and 
a v ~ 1.8 for the "M2" model. 

Too few data points exist up to now in or- 
der to constrain accurately many parameters from 
the models, such as the galaxy formation scenario 
l|Madau et alj fl998) (only good quality data at redshift 
> 1.5 could eventually distinguish between the two fa- 
vored galaxy formation scenar ios), the progenitor model 
( Rui z-Lapuente fc Canalll l998l). environ mental effects like 
the metallicitv faobavashi et alJl2000|) or the cosmology 
l|Dahlen fc Franssonlll999|) . 



9.3. Conclusion 



Fig. 15. Derivation of the evolution index a using all 
available SNIa rate measurements (see Table |3J), within 
the (0.3, 0.7) cosmological model. 



and the supernova explosion 6 , but also on the galaxy for- 



6 According to lDahlen fc Franssonl dl999l) . the r = 0.3 Gyr 
model corresponds to the double degenerate progenitor sce- 
nario, while the r = 1 Gyr model simulates the single degen- 
erate - or cataclysmic - progenitor scenario; the r = 3 Gyr 
expands the range to cover all likely models. 



We have measured the type la supernova explosion rate 
at a redshift of 0.13. Combined with other measurements 
at different redshifts, it suggests a slight evolution of the 
rate with the redshift. Nevertheless, the scatter and the 
still large error bars of the current measurements do not 
allow us to make definite conclusions about it. 

Future supernova searches will shed more light on the 
rate evolution. Among them, two dedicated supernova 
searches, are well suited to improve and measure the evolu- 
tion of the explosion r ate. The Nearby Supernova Factory 
l|Aldering et alJ 12002) is a nearby (z < 0.1) supernova 
search using two dedicated telescopes, one for the search 
and one for the follow-up of discovered events using an 
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Fig. 16. Expl osion rate evo l ution per unit of time and comoving volume. The set of curves labeled "Madau" stands for 
the models bv lMadau et al.l l )l998|) . computed within the hierarchical clustering galaxy formation scenario; r is the time 
interva l between the white dwarf formation and the supernova explosion (in Gyr). The curves labeled "Sadat" stand 
for the lSadat et al.l l|l998j) models, using two different star formation rate evolution models (Ml a nd M2, the latter 



accounts for a possible d ust ex tinction cor rection). From left to right are the mea sured values bv ICappellaro et alJ 
l|l999j) . lMadgwick et all l|2003jh this work. lHardin et all l|2000|) . iTonrv et all l|2003[) and lPain et alJ <|2002h . We plot 
also the redshift evolution, (1 + z) 175 , corresponding to the best fit using the six measured values. The effect of the 
evolution of the galaxy luminosity density {a.j B ) is shown as well since it is included in the rate per comoving volume 
unit evolution. All data points and model curves have been computed with the (0.3, 0.7) cosmological model, and 
/170 = 1. The normalisation of models is set by adjusting them to the data. 



integral field spectrometer. It will be very helpful to mea- 
sure the local supernova rate using a very homogeneous 
set of a few hundred events, and it will also help reducing 
the uncertainties on the supernova light curve distribu- 
tion. Another supernova search in progress uses the wide 
field Megacam camera on the CFHT 7 . Supernovae are 
discovered in a rolling search mode up to z < 0.9; it 
will be possible to use them for measuring the evolution 
of the explosion rate for type la supernovae at a 15 % 
level of statistical uncertainty (per bin of 0.1 in redshift). 
Concerning other supernova searches well suited to com- 
pute the explosion rate, we can quote the Sloan Digital 
Sky Survey performing a spectrosc opic detection of nearb y 
supernovae, as already mentioned l|Madewick et al . 2003); 
the Lick Observatory and Tenagra Observatory Supernova 



7 http://snls.in2p3.fr 



Searches 8 detecting nearby event using automatic tele- 
scopes; at intermediate redshifts (z ~ 0.4), the ESO 
Supernova Search 9 is designed to search for supernovae 
in order to derive the rate; in the distant Universe, the 
Great Observatories Origins Deep Survey Treasury pro- 
gram 10 is searching for very distant type la supernovae 
using the Hubble Space Telescope. 
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